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Abstract 

We prove an adiabatic theorem for the ground state of the Dicke model in 
a slowly rotating magnetic field and show that for weak electron-photon cou- 
pling, the adiabatic time scale is close to the time scale of the corresponding 
two level system-without the quantized radiation field. There is a correc- 
tion to this time scale which is the Lamb shift of the model. The photon 
field affects the rate of approach to the adiabatic limit through a logarithmic 
correction originating from an infrared singularity characteristic of QED. 
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I. INTRODUCTION 



In this work we investigate the relation between adiabatic theorems for models that, like 
QED, allow for the creation and annihilation of photons, and the corresponding quantum 
mechanical models where the electron is decoupled from the photon field. We study this 
problem in the context of a specific and essentially soluble model: The Dicke model fL3 



The corresponding quantum mechanical model is a two level system, such as a spin in an 
adiabatically rotating magnetic field, which is a basic paradigm of adiabatic theory M . 



In the usual quantum adiabatic theorem |I|[T8| the gap between eigenvalues plays an 



important role: It fixes the adiabatic time scale and determines the rate at which the 
adiabatic limit is approached. There is no such gap in the corresponding QED models so 
the nature of the adiabatic theorem in the two cases has qualitatively different features. For 
example, there is no gap in the spin-boson and Dicke models (for weak coupling) both of 
which describe a two level system in a radiation field. 

The first problem we address is whether there is an adiabatic theorem for the ground 
state in a radiation field. Assuming a positive answer, the second question is, what property 
of the QED model, plays the role of the gap in the adiabatic theorem. Another way of 
phrasing this question is how does the adiabatic time scale of the two level system compare 
to that of the QED model? Are the two close in the limit of small fine structure constant, 
a, and if so, how close? The third question compares the rate of approach to the adiabatic 
limit in the two models. 

Consider a two level system, such as a spin or a twofold Zeeman split atomic level, in 
an external magnetic field pointing in the z direction. When radiation effects are neglected, 
the corresponding Hamiltonian is 

H = ma z , m = fiB. (1) 
The corresponding Dicke model is 

H D = H ® 1 + a" 1 1 <g> E + ^faa + ® a\f) + ^fa cx_ ® a(/), (2) 

where 

E = f \k\a\k)a{k)d d k, (3) 



dQ 



an> 

/(*) = ^^i|{c-**,p}|^>, (4) 

with ipj{x) the atomic wave functions of the two level system. Polarization indices are 
omitted since the helicity of the photon does not play an interesting role in the questions we 



{•, •} stands for anticommutator 
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study. We use atomic units where e = h = lsoa = l/c = 1/137 is small, /i, the magnetic 
moment, is also of order a in these units. 

Following M. Berry's let us consider the case where the magnetic field changes its 
direction adiabatically and has fixed magnitude. The time dependent Hamiltonian for the 
two level system is 

H(s) = iiB(s) ■ a = U(s)HU*(s), (5) 

with U(s) G SU{2) the appropriate rotation. The corresponding adiabatic Dicke model has 
the time dependent Hamiltonian: 

H D {s) = (U(s) ® 1) H D (U*(s) ® 1) . (6) 

Our aim is to compare the adiabatic evolution of the ground state of H(s) with that of 
H D (s). 

Adiabatic theorems for quantum systems coupled to a field have been studied in p2|,|TT 



In |22[] Narnhofer and Thirring give characterization of extremal KMS states by adiabatic 
invariance. When applicable, this result says, in particular, that the ground state is adiabatic 
invariant. The characterization depends on a condition of asymptotic abelianess which 
does not hold for the models we consider. In JT1J Davis and Spohn give a derivation of 



linear response theory for a system coupled to a bath in the adiabatic limit. The notion 
of adiabaticity in this work is such that the coupling between the field and the quantum 
system vanishes in the adiabatic limit. This is not a standard notion of adiabaticity. 

Let us now describe our results. First, we show that there is an adiabatic theorem for 
the ground state of the Dicke model, even though the model has no spectral gap to protect 
the ground state. Second, we show that the the distance to a nearby resonance in the Dicke 
model plays the role of a gap. Third, we show that the adiabatic time scale for the Dicke 
model and the two level system agree in the limit of small a. The (inverse) of two time scales 
differ by the Lamb shift of the Dicke model. And finally, we show that the approach to the 
adiabatic limit in the two models is different: While in the two level system the approach to 
the adiabatic limit is with an error O (Jp), the approach to the limit in the Dicke model is 

with an error of O ( v/l ° g r ^ . The logarithm comes from an infrared divergence characteristic 



of QED. 

Although the results we derive here are for a rather special model we suggest that some- 
thing similar happens also for more realistic models. The success of the quantum in numerous 
applications that depend on a correct prediction of the adiabatic time scale is evidence that 
at least the time scale aspect of our results may well carry over to more realistic models. It 



would be interesting to know if this is indeed the case for the Spin-Boson model ||IT| , |r2| , |2T 
The spin-boson model is a more realistic QED version of a two level system, which, unlike 
the Dicke model, is not explicitly soluble. However, as much progress in the spectral analy- 
sis of the spin-boson problem has been recently made, the problem we pose here may be a 
reasonable challenge. 
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II. THE ADIABATIC THEOREM AND A COMMUTATOR EQUATION 



In this section we explain what we mean by "adiabatic theorem" , and give a condition 
for an adiabatic theorem to hold. This condition is that the commutator equation, Eq. (p~0|) 
below, has solutions X, Y which are bounded operators^. We also introduce notation, 
terminology, and collect known facts that we need. To simplify the presentation, we shall 
stay away from making optimal assertions. 

We consider Hamiltonians that are bounded from below, and choose the origin of the 
energy axis so that the spectrum begins at zero. Let H(s) > be a family of such self-adjoint 
Hamiltonians. The unitary evolution generated by the Hamiltonian, U T (s), is the solution 
of the initial value problem: 

i U T (s) = tH(s)U t (s), U t (0) = 1, s G [0, 1]. (7) 

t is the adiabatic time scale, and we are concerned with the limit of large r. The physical 
time is t = rs G [0, r). Since r is large H(s) = H(t/r) varies adiabatically. We assume that 
all operators are defined on some fixed dense domain in the Hilbert space. 

The (instantaneous) ground state is in the range of the kernel of H(s) and we assume 
that the kernel is smooth and one-dimensional. Let P(s) ^ be the projection on the 
kernel of H(s), i.e. H(s) P(s) = 0, dimP = Tr P = 1. By smoothness we mean that P(s) 
a bounded operator. 

The adiabatic theorems we consider are concerned with the large time behavior of the 
evolution of the ground state where t = 0(r) or, equivalentlys = 0(1). The smoothness 
of the kernel implies that there is a natural candidate for an adiabatic theorem for the 
ground state, which is independent of whether H(s) does or does not have a gap in it 
spectrum. Namely, that if ip(0) G Range P(0) at time s = 0, then it evolves in time so that, 
ip T (s) = U T (s) ip(0) lies in Range P(s) at time s in the adiabatic limit, r — > oo. 

To formulate the adiabatic theorem with error estimates we need to get hold of adiabatic 



phases 0. To do that we introduce the adiabatic evolution of Kato [fLSfl : Let Ua(s) be the 
solution of the evolution equation 

U A (s) = [P(s),P(s)]U A (s), U A (0) = 1, se[0,l]. (8) 

It is known that 

U A (s)P(0) = P(s)U A (s). (9) 

That is U A (s) maps Range P(0) onto Range P(s). We can now formulate the basic adia- 
batic theorem : 

Theorem II. 1 Let H(s)P(s) = for all < s < 1, with P differentiable projection 
on the ground state, with ||P(s)|| < D. Suppose that the commutator equation 



2 for X we also need that its derivative is bounded 
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[P(s),P(s)} = [H(s),X(s)}+Y(s), (10) 
has operator valued solutions, X(s) and Y(s) so that for e \ 

\\X{s)\\ + \\X{s)\\ < ^{|^| , \\Y(s)\\<Ce», (11) 

with /i, v > 0. Then 

||(^( 5 )-f/A( 5 ))P(0)||<C'| r 1 ^ r , sG [0,l]. (12) 

Remarks: 1. In the case that there is a gap in the spectrum, one can always find X(s) 
bounded so v = 0, and Y — 0, see 0. X, and therefor also C, is of the order of (gap) -1 . 
This gives error of 1/r, and generalizes the adiabatic theorem of Born and Fock and Kato 
for discrete spectra, to more complicated spectra provided there is a gap. 

2. The theorem says that the physical evolution clings to the instantaneous spectral 
subspace. In particular, if P is one dimensional, it says that the physical evolution of the 
ground state remains close to the instantaneous ground state. 

3. Here, and throughout, we are concerned only with the adiabatic theorem to lowest 
order. If s is chosen outside the support of P then much stronger results can be obtained. 
See e.g. |BJ . 



4. The adiabatic time scale To set by this theorem is tq = 0((2 + D)C). 
Proof: Let W(s) = U A (s)U T (s)), with W(0) = 1. From the equation of motion, and the 
commutator equation, Eq. (p!Q|). 

P(0)W{s) = -P(0)rf A (s)(iTH(s) + [P(s),P{s)])U T {s) 
= -Ui(s)P(s)(irH(s) + [P(s),P(s)])U T (s) 



= -U A (s)P(s) [P(s),P(s)]U T (s) 

= -U A (s) P(s) ([H(s), X(s)] + Y(a)) U T (s) 

= -U A (s) P(s) ( - X(s) H(s) + Y(s)) U T (s) 

= - P(0) U A (s) X(s) U T (s) - P(0) U\{s) Y{s) U T (s). (13) 
r 

To get rid of derivatives of U T , which are large by the equation of motion, we rewrite the 
first term on the rhs (up to the P(0) on the right) as : 



U\{s) X(s) U T (s) = (U\(s) X(s)U T (s)) - U\{s) X(s) U T (s) - U\{s) X{s) U T {s) (14) 
[u A (s) X(s) U T (s)) - U{(s) X(s) U T (s) + U\{s) [P(s), P(s)} X(s) U T (s). 
From this it follows, by integrating, that for s £ [0, 1] 

\\(U T (s) - U A (s)) P(0)|| = ||P(0)(f/t( S ) - U A (s)) || 

||P(0)(1 - W(s)) || < Ce» + (2 + r D) C | j ^ (15) 



Choosing e = r ^ gives 

\m(s)-U A (s))P(0)\\<cS [ T ^ • (16) 

This concludes the proof of the theorem. L~H 
It is convenient to rewrite this solvability condition in a way that one needs to solve for 
a fixed X and Y rather than functions X(s) and Y(s). This is accomplished by 

Corollary II. 1.1 Let P(s) be the family 

P(s) = V(s)PV\s), V(s) =exp(isa). (17) 
It is enough to solve for the commutator equation 

iK=[H,X]+Y, K = {a, P} — 2PaP, (18) 
for fixed X and Y so that for e \ 

wi £ c {id • » y » £ ^ < 19 > 

with /i, v > 0, and ||P(s)|| < X(s) and K(s) are then determined by the obvious unitary 
conjugation. 

Proof: Since P(s) = V(s) PV^(s), we have 

P(s)=tV(s)[a,P]V\s), (20) 

and 

[P(s),P(s)} = iV(s) [[a,P],P] V\s) 

= tV{s)({a,P}-2PaP)v\s). (21) 

□ 

III. AN ADIABATIC THEOREM FOR A THRESHOLD STATE: THE 

FRIEDRICHS MODEL 

As a warmup, and as a preparation for the analysis of the Dicke model, let us prove an 
adiabatic theorem for the Friedrichs model which has a bound state at the threshold of the 
continuum. 

There is an inherent difficulty in the situation of a bound state at threshold in general, 
and in the Friedrichs model |T^-[16| in particular, namely, that a bound state at threshold is 
not a stable situation. Under a small deformation of the Hamiltonian, the ground state will, 
generically, split away from the absolutely continuous spectrum and a gap develops. Since 
our aim is to study families related by a unitary, this problem does not appear. That is, we 
consider the family H F (s) = V(s) H F V\s) where H F has a bound state at threshold and 
V(s) is a smooth family of unitaries. 
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A. The Friedrichs Model 



We shall consider a family of Hamiltonians, closely related to the standard Friedrichs 
model parameterized by the scaled time s, a real number d > that plays the role of 
dimension, and a function / that describes the deformation of the family. Since we are only 
interested in the low energy behavior of the family we shall introduce an "ultraviolet cutoff" 
to avoid inessential difficulties. 

The Hilbert space of the Friedrichs model (with an ultraviolet cutoff) is TC = (D © 
L 2 ([0, 1], k d ~ x dk). A vector ip e H is normalized by 



0=1 ™i imi 2 = i/3| 2 + [inmv-uk, $ ec. (22) 



M w2 = 

We choose a special, and trivial, case of a diagonal Hamiltonian whose action on a vector ip 
is as follows: 

"*-(!!2)U,)-U,)- 

H has a ground state at zero energy with projection 

'-(is)- <*> 

The rest of the spectrum is the unit interval [0, 1], and is absolutely continuous. The density 
of states in this model is proportional to E d ~ x . 

We construct the family H(s) by conjugating H with a family of unitaries: 

V f (s) = expzscr(f), a(f) = ( ° /} { Q fl ^ , (25) 

where / is a vector in L 2 ([0, 1], k d ~ x dk). 



Theorem III.l Let Hp(s;d, f) be the family of Friedrichs models with a ground state 
at threshold for all s 

H F (s;d,f) = V f (s)H F V}(s). (26) 

Suppose that 

g(k) = i AT 1 /(*) G L 2 ([0, 1], k d - x dk\ V f (s) =expts a(f) , (27) 

then the quantum evolution of the ground state of H F (s; d, f) is adiabatic and its deviation 
from the instantaneous ground state is, at most, 0(l/r). 

Remarks: 1. Note that if the conditions in the theorem hold in dimension do, then they 
hold in all dimensions d > do. The physical interpretation of that is that the density of 
states at low energies decreases with d. So, even though there is spectrum near zero, there 
is only very little of it. 



7 



2. If is not in L 2 there may still be an adiabatic theorem with slower falloff in r by 
accommodating 7^0. An example will be discussed in the next section. 

3. The Friedrichs model is vanilla: Hp has no interesting energy scale to fix the adiabatic 
time scale. The scale is set by the perturbation alone: T = O((l + ||/|| 2 )||g||). This is quite 
unlike the case in the usual adiabatic theorem and unlike what we shall show for the Dicke 
model. 

Proof: In this case K of Corollary 2.1 is K = cr(f). With g G L 2 , a(g) is a bounded (in 
fact, finite rank) operator and an easy calculation gives 

[H F ,a(g)]=(? kg) { ~ k9 l)=ia(f). (28) 

Hence 

X = a(g), Y = 0, (29) 

solve the commutator equation, Eq. (|18D, with a bounded X(s) and Y(s) = 0. d 



IV. ADIABATIC THEOREM FOR THE DICKE MODEL 



In this section we describe an adiabatic theorem for the Dicke model [I3| that says that 
the an adiabatic rotation of a two level system evolves the ground state so that it adheres to 
the instantaneous ground state and the time scale, at least in three dimensions, is essentially 
the time scale fixed by Quantum Mechanics without photons. The rate of approach to the 
adiabatic limit is different from that of a two level system and has a logarithmic correction 
in three dimensions. This section also collects known facts about the Dicke model that we 
need. 



A. The Dicke Model 



The Spin-Boson Hamiltonian in the canonical QED version of a two level system 
P~T| , [T2|J2T[] . The Dicke model is a simplified version of the Spin-Boson Hamiltonian in the 



rotating wave approximation. The rotating wave approximation, can indeed be motivated 
in the single-mode Dicke model. In the multi-mode case we consider the rotating wave ap- 
proximation is a name that describes which terms in the Spin Boson Hamiltonian are kept 
and which are not. 

The model describes a two level system coupled to a massless boson field in d dimensions. 
The Hamiltonian is: 

H D (m, d, f, a) = m (1 - P) <g> 1 + a' 1 1 g> E + <Ja a + <g> a f (/) + y/a a_ <g> a(f), (30) 

acting on the Hilbert space C 2 <S> T with T being the symmetric Fock space over L 2 (R d , d d k). 
Here 



P = nn . °+= nn > ° = in E = \k\a\k)a{k)d d k. (31) 
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m > is the gap in the quantum Hamiltonian (without photons). a(f) and a'(f) are the 
usual creation and annihilation operators on T obeying the canonical commutation relations 



[«(/), at (<?)] = </!<?> ■ 



(32) 



We denote by |0) the field vacuum and by Q the projection on the vacuum. 

It may be worthwhile to explain where the various powers of a in if come from. For the 
radiation field the a -1 comes from hu = hc\k\ which explains why the field energy comes 
with a large coupling constant. The yfa. has one inverse power of c from minimal coupling, 



e 

2ni<: 



(p-A + A-p). Half a power of a comes from the standard formula for the vector potential 

(33) 



A(x) :-- 



d A k 



I27TC 



\k\ 



(e- ik - x a\k) + e ik - x a{k)). 



Compare e.g. p0 |. 

With reasonable atomic eigenfunctions, f(k), Eq. (|j) has fast decay at infinity and the 
model is ultraviolet regular. In the infrared limit f(k) behaves like 



/(*)-> -M/?U / (^)(W> 2 )(x) -{V^ 1 )*{x)^ 2 {x))d d x. 



In particular we see that for small k 



/(*) = K 




(34) 



(35) 



The square root singularity is a characteristic infrared divergence of QED, and it has con- 
sequences for the adiabatic theorem as we shall see. Note that with / having a square root 
singularity the model makes sense (as an operator) provide d > 1, for otherwise a/(f) is ill 
defined since / is not in L 2 . 

An important parameter in the model is 



£ = {f 



l 



k 



f 



Bearing in mind the square root singularity of / we see that 

, d d k 



\k 



2 : 



(36) 



(37) 



is finite for all d > 2. 



B. Spectral Properties 



What makes the Dicke model simple is that it has a constant of motion [0. If we let 
N = J a){k) a{k)d d k be the photon number operator, then M commutes with H D where 

N= ° N+1 ^j=l®N + P®l. (38) 

The spectrum of Af is the non-negative integers. The spectral properties of H D (m, d, f, a) 
can be studied by restricting to subspaces of Af. 
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a. Af = : The kernel of TV is one dimensional and is associate with the projection 




p= „ ■ (39) 



Q is the projection on the field vacuum. It is easy to see that P Hc(m, d, f, a)P = 0, so the 
model always has a state at zero energy. This state may or may not be the ground state. It 
is the ground state if a 2 £ < m |17| . 



b. Af — 1 : The space is basically 7i of the Friedrichs model. The correspondence of 
vectors in the two spaces is 

(40) 





The Hamiltonian action in the Friedrichs model language is: 

ff »M/,«)»(^|f fl ). (41) 

It is a standard fact about the Friedrichs model [0,0] that provided 

a 2 S < m, (42) 

the model has no bound state, and the spectrum is [0, oo) and is absolutely continuous. Since 
/ has square root singularity at the origin, (and has fast decay at infinity), this condition 
holds for d > 3 if a (or /) is not too large. In three dimensions, provided that the level 
spacing m » a 2 in atomic units, (about 10 -3 eV), the inequality holds. In two dimensions 
the left hand side is log divergent, and the spectrum in the Af = 1 sector has a bound state 
at negative energy. This state lies below the bound state of the Af = sector. We do not 
consider this situation and henceforth stick to d > 3. 

c. Af > 2 : It is known [0] that the bottom of the spectrum in all these sectors is at 



zero if (E2) holds. 



C. Adiabatic Rotations 

Suppose, that the two level system of the Dicke model describes e.g. two Zeeman split 
energy levels of an atom in constant external magnetic field B pointing in the z direction. 
Rotations about the z axis do not change the orientation of the magnetic field, and commute 
with Af and are uninteresting. Rotations about the x axis change the orientation of the 
magnetic field and are implemented by 

V(s) = exp(isa x ) 1. (43) 

Such rotations do not commute with Af. Indeed, 

[Af,a] = ( _° J ® 1 = J® 1, a = o x ®l. (44) 
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D. The Adiabatic Theorem 



Theorem IV. 1 Let H D (s;m, d, /, a) = V(s)H D (m, d, /, a)V"t(s), s G [0,1] be the family 
of time dependent Dicke models with / square integrable, with square root singularity at 
k = 0; m > a 2 (f\ m |/); d > 3 and V(s) = exp (i sa) as in Eq. (|43"D . Then, C/^, the adiabatic 
evolution associated with the ground state of H D (s;m,d,f,a), and U T , the Schrodinger 
evolution are close in the sense that 



\\(U A (s)-U T (s))P(0)\<C\ 



if d > 3 
d = 3. 



(45) 



The time scale is determined by m — a 2 £ and coincides with the gap without photons, m, 
up to a correction by the Lamb shift, a 2 S. 

Proof: From Corollary 2.1 we find K = a ® Q. We will first show that a solution of the 
commutator equation, Eq. (|i~8"D, for ci > 0, is 



where 



x = ^M r = o, 

m — cr 6 



X 1 = J®fi, X 2 (</) = (a f (#)ft + h.c), 5 



l a 2 777 



(46) 



(47) 



Note that the gap of the two level system m is renormalized to m + i a (f\g), which is just 
the Lamb shift (See appendix). This is a small correction, of order a 2 . 
A useful formula we shall need is 



Ea!(g)n = a ] {\k\g)VL. 
Let us compute the commutators of Xi, X 2 with H: 



(48) 



[H,X X ] 



E 
a 



a ofi {f) \ 
aa(f) m+f ) 



-n 
n o 



(49) 



For the second commutator 

E 



[H,X 2 



_ ao f (/) 
aa(f) m + f 



a^{g) Q + h.c 




a 



P ® {a){\k\g)Q, - Qa(\k\g) + ^ 





(f\9) 



-{g\f)\ 



n. 



So, if we take g of Eq. fl47|) then 

i Z x - X 2 ] = i (m - a 2 £) a ® a 



(50) 



(51) 



We see that we can formally solve for the commutator equation, Eq. (|I8]) provided £ is finite. 
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This is, however, not the only condition. X is a bounded operator in the Hilbert space 
provided g G L 2 , for otherwise a'(g) is ill defined: 



rX^d d k ~ / -^d d k < oo. (52) 

The integral is finite if d > 4 but is logarithmically divergent in d = 3. For d = 3 we need 
to squeeze X back to the bounded operators. We do that by allowing for 7^0. 

Let Xe be the characteristic function of a ball of radius e and x% = 1 — Xe and let g c e = \%9 



and g £ = Xe9- Let us take X 2 (g £ ), which is well defined and its norm is 0{a^ J\ log e | ) . For 



X we take, as before, 



x= lX -:^'' ] , 

m + i^(f\gf) 



From this 



_ 3 
1 + «2 


\oge 


1/2 




m - 


-a 2 £ 





X\\=0\ ' 1 » ■ (54) 



For K we take 



(m + tV^(f\9e))Y=[H,X 2 (g)-X 2 (g c £ )] = [H,X 2 (g £ )] 

= - P ® (aU\k\g e )tt - tta(\k\g 
a v 




and we used the computation of the commutator Eq. (|50D . With / having a square root 
singularity, 

y |m — a z c | / 

This puts us in the frame of theorem II. 1, except for the minor modification the log appears 
with a square root. Chasing the square root establishes the main result. L~H 
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APPENDIX A: RESONANCE AND LAMB SHIFT OF THE DICKE MODEL 

The M = 1 sector of the Dicke model has a resonance that serves to define the Lamb shift. 
The resonance is a solutions of the analytically extended eigenvalue equation, Eq. (|A1|), see 



I4] , p!6|pl| . The real part of the shift is, by definition, the Lamb shift of the model, while the 
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imaginary shift is the life time. For d > 3, the Lamb shift is dominant and the life time is a 
higher order in a. For the application to the adiabatic theorem we need only the dominant 
contribution, i.e. only the Lamb shift. Computing the Lamb shift is easy. Computing the 
life time is harder. For the sake of completeness we compute both, even though we only 
need one. 

The eigenvalue equation is 

E — m = a 2 G(aE), (Al) 
where G(e) is defined as the analytic continuation from the upper half plane of 

G(e)= [ -^-d d k, 3e>0. (A2) 
Jn d e — \k\ 

By taking the imaginary part, it is easy to see that Eq (|A1|) has no solution in the upper 
half plane. To solve the equation in the lower half plane one needs an explicit expression, at 
least for small a, and e near am, of this analytic continuation. Then, we can solve Eq. (|A1|) 
by iteration, and to lowest order we have 

E r « m + a 2 G(am). (A3) 

Clearly G(am) — > —£, in the limit a — > 0, so to leading order 

E r ^m-a 2 £. (A4) 

To this order, one does not see the imaginary part of the resonance energy, a 2 £ is, by 



definition, pOfl , the Lamb shift of the model. It may be worthwhile to point out that the 
Lamb shift for the Hydrogen atom, ||, is actually of higher order, namely, a 3 log(« _1 ). Since 
the Lamb shift of Hydrogen also involves an ultraviolet regularization, while the present 
model is ultraviolet regular, it is not surprising that the order of the two is different. What 
is surprising is that the order of Hydrogen is higher rather than lower. 

Estimating the life time is, as we noted, irrelevant to the adiabatic theorem. So a reader 
will loose little by skipping the rest of this Appendix. However, for the benefit of the reader 
who is interested in how the computation of the life time goes, it is given below. 

We shall show below that for d > 3, and |e — am\ < am, the analytic continuation of 
G(e) to the lower half plane, and to the next relevant order in a, is given by 

G(e) = — £ — in KQ d e d - 2 , 3 e < 0, (A5) 



where K is as in Eq. fl3"5]), and Cl is the surface area of the unit ball in d dimensions. From 
Eq. (|A3|), and taking into account Eq. (1351), we get for the Lamb shift and the life-time: 



E r « m - a 2 £ - ia 2 n K Vt\ma) d - 2 

= m-a 2 £- ia 2 n£l d (ma) d - 1 |/(am)| 2 . (A6) 

The life time is higher order in a than the Lamb shift, and is of order a d . For d = 3 this is, 
indeed, the order of the life time of atomic levels that decay by dipole transition. For small 
a the Lamb shift dominates the life time, both in the Dicke model and in Hydrogen. 
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It remains to show that the analytic continuation of G(e) to the lower half plane in a 
neighborhood of ma, is indeed given by Eq. ( |A5|) . This can be done as follows: Let B r be 
a ball of radius r = 2ma. Then, in the upper half plane 

G(e)=([ + [ )^-d d k = G r (e) + G c r (e). (A7) 
\JBr JB$) e — \k\ 

Clearly, G c r {e) extends analytically to a half circle in the lower half plane |e — am\ < am. 
In the limit of a — > 0, by continuity, 

G c r (0) -> -£. (A8) 

This is the dominant piece, and it is real. 

Consider the analytic continuation of G r (e) for \e— am\ < am. Since, for small argument 
f(k) is given by Eq. (p5|), one has (in the upper half plane) 

GJe) = Ktt d / dk = KQ d / dk, (A9) 

jo e — k J<y e — k 

where 7 is the obvious semi-circle in the complex k plane and Q d the surface area of the unit 
ball in d dimensions. The right hand side is analytic in e in the lower half plane provided 
|e — am I < am, and so gives the requisite analytic continuation. Since e is small, and of 
order a, to leading order, we have 

GJe) = KQ d I (fc ~ e+ , e)d ' 2 dk 
J~i e — k 



= -KQ d J2 y j J e^" 2 J{k - e)'- 1 dk 
w -KQ d e d - 2 I 



J 

dk 



k — e 



KVt d e 



d d-2 



(1am\ \ 
in +logf— —j + 0(aloga)J (A10) 



and the error term in approximation that we did not compute is real and being sub-dominant 
to 8 is irrelevant. 
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